# Exercicio 7 - Regresses Lineares Simples
# Camila Galheigo Coelho 


setwd("C:/Users/Camila/Documents/PESQUISA/1_DOUTORADO/Disciplinas/R")
getwd()

## Altura na infncia e na vida adulta ##

h.adu<-c(71,63,63,67,68,68,70,64)
h.inf<-c(39,30,32,34,35,36,36,30)
#Q1.
mod.h.fxet<-lm(h.adu~h.inf)
mod.h.fxet
summary(mod.h.fxet)
anova(mod.h.fxet)

# existe uma relao significativa entre a idade aos dois anos de idade e na idade adulta,
# conforme confirma o p = 0,00386

#Q2.
plot(h.adu~h.inf,xlab="Altura adulto", ylab="Altura infante")
abline(mod.h.fxet, col="blue")

h.adu.crenca<-2*h.inf
mod.h.fxet.cre<-lm(h.adu.crenca~h.inf)
summary(mod.h.fxet.cre)
anova(mod.h.fxet.cre)
abline(mod.h.fxet.cre, col="red")
# as linhas de regresso dos modelos apresentam inclinaes distintas e portanto preveem relaes distintas entre a altura na infncia e na vida adulta.   

#Q3. 
confint(mod.h.fxet) #calcula os intervalos de confianca dos coeficientes
coef(mod.h.fxet)

#Q4. 
confint(mod.h.fxet.cre) #os valores no esto dentro dos intervalos de confianca da regresso linear simples.
coef(mod.h.fxet.cre)



## Seriemas e Carcars ##
aves.cerrado<-read.table("aves_cerrado.csv",header=TRUE,sep=";")
str(aves.cerrado)
head(aves.cerrado)
aves.cerrado$carcara[is.na(aves.cerrado$carcara)]<-0  #substituindo NA por zero
aves.cerrado$seriema[is.na(aves.cerrado$seriema)]<-0  
aves.cerrado$urubu[is.na(aves.cerrado$urubu)]<-0
table(aves.cerrado$fisionomia)
aves.cerrado$fisionomia[aves.cerrado$fisionomia=="ce"]<-"Ce"  #arrumando o Ce que estava em minsculo
table(aves.cerrado$fisionomia)
aves.cerrado$fisionomia <- factor(aves.cerrado$fisionomia, levels=c("CL","CC","Ce"))
str(aves.cerrado)
summary(aves.cerrado)
head(aves.cerrado)

#Q1.

par(mfrow=c(1,3))
# CL - campo limpo
plot(aves.cerrado$seriema[aves.cerrado$fisionomia=="CL"]~aves.cerrado$carcara[aves.cerrado$fisionomia=="CL"], xlab="caracara", ylab="seriema", main="Campo limpo")
mod.CL<-lm(aves.cerrado$seriema[aves.cerrado$fisionomia=="CL"]~aves.cerrado$carcara[aves.cerrado$fisionomia=="CL"])
abline(mod.CL)

##Ce - cerrado
plot(aves.cerrado$seriema[aves.cerrado$fisionomia=="Ce"]~aves.cerrado$carcara[aves.cerrado$fisionomia=="Ce"], xlab="caracara", ylab="seriema", main="Cerrado")
mod.Ce<-lm(aves.cerrado$seriema[aves.cerrado$fisionomia=="Ce"]~aves.cerrado$carcara[aves.cerrado$fisionomia=="Ce"])
abline(mod.Ce)

##CC - campo cerrado
plot(aves.cerrado$seriema[aves.cerrado$fisionomia=="CC"]~aves.cerrado$carcara[aves.cerrado$fisionomia=="CC"], xlab="caracara", ylab="seriema", main="Campo cerrado")
mod.CC<-lm(aves.cerrado$seriema[aves.cerrado$fisionomia=="CC"]~aves.cerrado$carcara[aves.cerrado$fisionomia=="CC"])
abline(mod.CC)

#aparentemente no Campo Limpo quanto mais Caracars so avistados menos Seriemas o so.

#Q2.

summary(mod.CL) #p=0,105
summary(mod.Ce) #p=0,245
summary(mod.CC) #p=0,645
#no h relao entre nenhuma das fisionamias dado que nenhuma das regressoes  significativa.

## Resduo Iris  ##
library(datasets)
iris
setosa<-iris[iris$Species=="setosa",]
setosa #objeto "setosa" contendo apenas dados referentes  especie Iris setosa

#Q1.
larg.comp.sepala<-lm(setosa$Sepal.Width~setosa$Sepal.Length) 
larg.comp.sepala
plot(setosa$Sepal.Width~setosa$Sepal.Length,xlab="Comprimento spalas (cm)",ylab="Largura spalas (cm)",main="Iris setosa")
abline(larg.comp.sepala)
anova(larg.comp.sepala)

#Q2. 
res1<-lm(setosa$Sepal.Width~setosa$Petal.Length)
res1
plot(setosa$Sepal.Width~setosa$Petal.Length,xlab="Comprimento petalas (cm)", main="Iris setosa", ylab="Largura das sepalas (cm)")
abline(res1)
residuos1<- residuals(res1)

res2<-lm(setosa$Sepal.Length~setosa$Petal.Length)
res2
plot(setosa$Sepal.Length~setosa$Petal.Length,xlab="Comprimento ptalas (cm)",ylab="Comprimento spalas(cm)", main="Iris setosa") 
abline(res2)
residuos2<- residuals(res2)

desc.efeito<-lm(residuos1~residuos2)
plot(residuos1~residuos2,xlab="Resduos do modelo Comprimento spala ~ Comprimento ptala",ylab="Resduos do modelo Largura spala ~ Comprimento ptala")
abline(desc.efeito)
anova(desc.efeito) 
 
#a regresso dada pela relacao entre as variaveis continua significativa (p=1.431*10^-09), apesar da correo para o comprimento. 

save.image()